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epochs between 321 and 1799 days after the explosion show three distinct ve- 



locity components. The narrow lines come from circumstellar gas and show 
both low and high ionization. This component has a low filling factor, and is 
photoionized by X-rays from the shock. The intermediate component, which 
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2500-5000 km s , and most likely comes from the ejecta. The hydrogen- and 
helium-dominated gas has a low ionization, a high density, and velocities that 
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extend out to >10,000 km s" 1 . Strong signatures of Lya-pumped fluorescence 
lines of Fe II are seen in the near-infrared and ultraviolet. The He/H ratio, ~ 0.3 
by number, and the nitrogen overabundance provide strong evidence for CNO 
burning products. The fluxes of the broad hydrogen and helium lines decrease 
considerably faster than the oxygen lines. The Ha line profile shows strong evo- 
lution, with the red wing decreasing faster than the blue. Possible scenarios, 
involving either a clumpy circumstellar medium, or an aspherical distribution of 
the surrounding gas, are discussed based on the line profiles and physical con- 
ditions. Finally, we propose that Type Iln supernovae have their origin in red 
supergiants in a superwind phase. 

Subject headings: stars: circumstellar matter — stars: mass loss — supernovae: 
individual: SN 1995N 



1. INTRODUCTION 

SN 1995N in MCG-02-38-017 (also known as Arp 261) was discovered by Pollas (1995) 
on 5 May 1995 (UT dates are used throughout this paper). Its approximate distance is 
24 Mpc, based on H = 75 km s" 1 Mpc" 1 and the measured recession velocity of the 
host galaxy (1834 km s -1 ; Theureau et al. 1998). Optical spectra of the object obtained 
by Benetti, Bouchet, & Schwarz (1995) and by Garnavich & Challis (1995) showed it to 
be a peculiar Type II supernova, with Ha emission (full-width at half-maximum intensity 
[FWHM] « 1500 km s" 1 ) and narrow lines of He I, He II, [N II], [O I], [O III], [Ne III], 
[Ca V], Fe II, [Fe II], [Fe III], [Fe VII], and [Fe X]. The intensity ratios of the narrow oxygen 
lines suggest a high electron density, > 10 6 cm" 3 (Garnavich & Challis 1995). Benetti et 
al. (1995) estimated it to be at least 10 months old, so we will assume an explosion date of 
4 July 1994, befitting for a supernova (in the opinion of the U.S. coauthors of this paper). 
Radio emission from the supernova was discovered at 3.6 cm on 16 June 1995 (Van Dyk et 
al. 1996). X-ray emission was discovered with ROSAT on 23 July 1996, and the source was 
subsequently observed with ASCA on 20 January 1998 (Fox et al. 2000). The unabsorbed 
X-ray luminosity was 1 x 10 41 erg s -1 in the 0.1-10 keV band. 

The narrow emission lines observed from SN 1995N place it in the Type Iln category 
(Schlegel 1990; Filippenko 1991a,b, 1997), which often show a slow optical decline compared 
to other Type II supernovae (SNe II). Baird et al. (1998) have noted the slow evolution of SN 
1995N. The high X-ray and radio luminosities of SN 1995N place it in a subset of SNe Iln 
whose other well-studied members are SN 1986J (Rupen et al. 1987; Weiler, Panagia, & 
Sramek 1990; Leibundgut et al. 1991; Houck et al. 1998), SN 1988Z (Stathakis & Sadler 
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1991; Filippenko 1991a,b; Turatto et al. 1993; Van Dyk et al. 1993; Fabian & Terlevich 
1996; Aretxaga et al. 1999), SN 1978K (Ryder et al. 1993; Chugai, Danziger, & Delia Valle 
1995; Schlegel, Petre, & Colbert 1996), and SN 1998S (Leonard et al. 2000, Gerardy et al. 
2000; Fassia et al. 2001). All of them have X-ray luminosities >10 40 erg s -1 , which are at 
least an order of magnitude higher than the X-ray luminosities of other SNe II. These are 
also among the most luminous radio SNe, although SN 1979C does rival them in luminosity, 
as does SN 1998bw, a very different kind of object and at a much earlier epoch (e.g., Galama 
et al. 1998; Kulkarni et al. 1998). 

The early optical observations of SN 1988Z (Stathakis & Sadler 1991) showed velocities 
up to ~ 20, 000 km s _1 , and very long baseline interferometry (VLBI) radio observations of 
SN 1986J (Bartel, Shapiro, & Rupen 1989; Bartel et al. 1991) implied expansion velocities 
of ~ 15,000 km s _1 . At the same time, these SNe showed Ha emission lines with FWHM 
« 500 - 2000 km s~\ and at least SN 1988Z exhibited very narrow emission lines whose 
intensity ratios indicate high gas densities. This combination suggests that the rapidly 
moving supernova material is running into dense circumstellar clumps and is driving slower 
shock fronts into them (Chugai 1993; Chugai & Danziger 1994; Chugai et al. 1995). In 
addition, the combination of high velocities with strong circumstellar interaction requires 
that the supernova have considerable high-velocity gas in order not to be rapidly decelerated. 
For SN 1988Z, Chugai & Danziger (1994) estimate that the total mass in supernova ejecta 
is only ~ 0.5 M if the explosion energy is the standard 10 51 ergs. 

In order to more fully understand the nature of this peculiar subset of SNe Iln, we 
have undertaken a detailed spectroscopic study of SN 1995N, one of the most important 
and thoroughly studied examples. We present the ground-based and Hubble Space Telescope 
(HST) observations in § 2. In § 3, we discuss the three main kinematic components visible in 
the emission-line profiles. Synthetic spectra are presented in § 4, and the temporal evolution 
of features is examined in § 5. We discuss the results and their implications in § 6, and we 
summarize our conclusions in § 7. 

2. OBSERVATIONS 

Figure 1 shows a 3-minute i?-band image of MCG-02-38-107 and SN 1995N taken with 
the FORS1 instrument on the ESO Very Large Telescope (VLT) on 11 May 1999 (UT dates 
are used throughout this paper). Figure 2 is an HST WFPC2 /-band (F814W) image of SN 
1995N obtained on 22 July 2000; the supernova is easily visible at ~ 20.5 mag. These two 
images demonstrate that the supernova occurred in a region relatively free of background 
emission and dust. 
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The main observations of SN 1995N presented in this paper consist of ground-based 
optical spectra obtained on nine epochs and ultraviolet (UV) HST spectra on one epoch. 
The log of the observations is given in Table 1, where the wavelength ranges, UT dates, and 
approximate epoch since explosion are summarized. The assumed explosion date of 4 July 
1994 is uncertain, but this does not appreciably affect our results because our data were taken 
long after the explosion; the spectra cover epochs of 321-1799 days. Standard procedures 
were followed for the spectral extractions and reductions, which included sky subtraction, 
wavelength and flux calibration, and removal of telluric absorption (e.g., Matheson et al. 
2000). 

The ground-based observations fall in two main groups: within ~ 2 years after dis- 
covery, and ~ 4 years after discovery. The former were obtained primarily with the Kast 
spectrograph (Miller & Stone 1993; 2" slit) on the 3-m Shane reflector at Lick Observatory, 
but the Low Resolution Imaging Spectrometer (LRIS; Oke et al. 1995; 1" slit) on the Keck-I 
telescope was used to acquire an excellent spectrum (day 1007) nearly contemporaneously 
with one of the Lick spectra (day 1012). All of the Lick and Keck spectra (except day 1012) 
are shown in Figure 3. One can see that the spectrum evolved slowly during the time interval 
shown. In this paper, we present a detailed analysis of only two of these spectra (days 716 
and 1007, which are among the two best ones). 

The very late-time (~ 4 years) ground-based spectra were taken through a 1" slit with 
FORS1 on the VLT at ESO. The day 1769 spectrum was obtained with gratings 300V and 
3001, giving a resolution of ~ 13 A. For the day 1799 spectrum the 600B, 600R, and 6001 
grisms were used, giving a higher resolution of ~ 5 A. The day 1769 observing was done 
in good conditions, while the data from day 1799 were taken through clouds. The latter 
spectrum has higher signal-to-noise ratio (S/N) than the former, however, so emission-line 
measurements are reported for it (scaled to the flux level of the day 1769 spectrum). 

The HST spectra on day 943 were obtained with the Faint Object Spectrograph (FOS) 
with the G400H, G270H, and G160L gratings. Total exposures times were 2860, 4770, and 
2110 s, and the resolution was 3.0, 2.1, and 6.9 A per diode, respectively. The spectra were 
calibrated by the FOS pipeline at the Space Telescope Science Institute. 

All spectra, except the VLT and HST spectra, were calibrated to the y-band and R- 
band photometry of Baird et al. (1998). The Keck spectrum (day 1007) has the highest 
S/N, and in Figure 4 we show it together with the HST spectra (day 943). The Keck and 
HST spectra were taken at different epochs, but from the day 1037 Lick spectrum we find 
a smooth evolution, and we normalize the flux level of the overlapping region between the 
HST G400H spectrum and the Keck spectrum to that of the Keck epoch. 
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The Dirbe/IRAS maps (Schlegel, Finkbeiner, & Davis 1998) give E B -v = 0.11 mag for 
the Galactic reddening of SN 1995N. In addition, there could be internal reddening from 
both the host galaxy and the circumstellar environment. The host galaxy is a blue compact 
galaxy, a class which is generally dust-poor, although H II regions can show substantial 
reddening (Kong & Cheng 1999). The fact that the supernova is located in the outskirts of 
the galaxy (Figures 1, 2) also makes it natural to expect a relatively low extinction. From 
the relative strengths of the higher members of the H I and He I recombination series, there 
is no indication of strong reddening (§ 4.1). The normal intensity ratio of UV O III] lines 
relative to the optical also indicates this. In the following, we use E B -v = 0.11 mag for the 
total reddening. 



3. KINEMATIC COMPONENTS 

Figure 4 shows that there are lines spanning a large range of velocities, and we can 
distinguish at least three kinematic components in the supernova. The narrow, unresolved 
lines are clearly distinct, while the other two components are less obvious and are shown in 
the Mg II AA2796, 2803, [O III] AA4959, 5007, and Ha lines at 943-1007 days on the same 
velocity scale in Figure 5. The Ha and [O III] profiles are from the Keck day 1007 spectrum 
and the Mg II profile from the HST spectrum. The Ha and [O III] lines have quite different 
profiles, with the former having a fairly peaked profile, but with weak wings reaching very 
high velocity, ~ 10, 000 km s _1 , while the latter has a more boxy profile reaching a velocity 
of ~ 5000 km s _1 . Another strong indication of two separate high- velocity components is 
the different temporal evolution of the [O I], [O II], and [O III] lines compared to the others 
(§ 5). Here we discuss these three components individually. 



3.1. Narrow Lines 

The narrowest lines in the spectra are unresolved, implying FWHM <500 km s _1 . 
They are listed in Table 2, where it can be seen that the emission is dominated by ions with 
moderately high ionization. Some of the lines provide useful diagnostics for the emitting 
region, and we begin by considering them. 

[O III]. The [O III] AA4959, 5007 lines are by far the strongest narrow lines and, 
together with [O III] A4363, provide a useful constraint on the electron temperature T and 
the electron density n e . The relatively low ratio of /(A4959 + A5007)/7(A4363) = 5.4 implies 
that collisional deexcitation of the 1 D 2 level is important, so that the ratio depends on 
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density as well as temperature (e.g., Filippenko & Halpern 1984; Osterbrock 1989). The 
densities corresponding to various assumed temperatures are in the range 7.3 x 10 5 < n e < 
7.5 x 10 7 cm" 3 for 0.8 < T 4 < 4 (or 7.3 x 10 5 < n e < 2.0 x 10 6 cm" 3 for 0.8 < T 4 < 2), 
where T 4 = T/(10 4 K). On day 1799 the VLT spectrum shows a similar value, /(A4959 + 
A5007)//(A4363) = 6.5, which gives 7 x 10 5 < n e < 4 x 10 6 cm" 3 for 1 < T 4 < 2. 

O III Bowen lines. In the near-UV spectral region, several Bowen fluorescence lines are 
seen, most notably O III AA3047, 3133, 3340, 3444. We measure J(A3047)//(A3132) « 0.30 
and J(A3444)//(A3132) « 0.25. Based on the Bowen process, Kallman & McCray (1980) 
give ratios of 0.18 and 0.28 for these lines, respectively. Considering that the uncertainties 
in the fluxes of both lines are at least a factor of two, the agreement is satisfying. The A3340 
line is severely blended with [Ne III] A3342 and [Ne V] A3345, and no reliable flux estimate 
can be obtained. 

The ratio 1(0 III A3132)//(He II A4686) can be used to estimate the Bowen yield for 
He II A304. Equation (25) in Kallman & McCray (1980), together with 7(A3444)/7(A3132) = 
0.28, gives a Bowen yield of y B = 0.28 x 1(0 III A3132)//(He II A4686). Correcting for 
reddening, we find 1(0 III A3132)//(He II A4686) = 1.05, and therefore y B « 0.3 for SN 
1995N at day 1007. As a comparison, Schachter, Filippenko, & Kahn (1990) find a mean 
yB ~ 0.24 for a sample of Seyfert nuclei, while Kallmann & McCray find that planetary 
nebulae have y B ~ 0.37. The dispersion in both values is considerable. 

[O II]. The spectra show that [O II] AA3726, 3729 are not detected, but that [O II] 
AA7320, 7330 are weakly present with narrow components. The fact that J(A3726+A3729)//(A7320+ 
A7330) < 1 implies that n e > 3 x 10 4 cm" 3 , provided that T < 30, 000 K. 

[N II]. The auroral [N II] A5755 line is present, but the nebular lines are superposed 
on the very strong Ha emission and it is not possible to obtain strong limits on them. 
We estimate /(A6548 + A6583)/7(A5755) < 12, which implies that n e > 6 x 10 4 cm" 3 for 
T « 20,000 K (Osterbrock 1989). 

[Fe VII]. The observed lines of [Fe VII] are AA3586, 3759, 5159, 5276, 5721, 6087. As 
discussed by Nussbaumer & Storey (1982) and Keenan & Norrington (1987, 1991), several of 
the [Fe VII] lines provide diagnostic information on the emitting region. There is temperature 
sensitivity of the lines, and collisional deexcitation causes transitions from the low-density 
limits to the high-density limits in the range n e pa 10 6 — 10 8 cm~ 3 . 

The observed, reddening-corrected ratio J(A3759)//(A6087) w 0.83 in the Keck spec- 
trum (day 1007), and ~ 1.1 in the 12 June 1999 VLT spectrum (day 1799). From Figure 2 
of Keenan & Norrington (1987), this sets a limit T > 17, 000 K for any density on day 1007, 
and T > 22, 000 K for day 1799. 
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In the day 1799 spectrum, with the highest resolution, there are three distinct peaks 
at the positions of [Fe VI] A5146, [Fe VII] A5159, and [Fe VI] A5176. Because of blending 
of these lines, it is difficult to obtain accurate fluxes, and depending on where we estimate 
the base of the line we get a ratio in the range 0.20 < J(A5159)//(A6087) < 0.73 and 
0.25 < J(A5159)//(A6087) < 1.1 on days 1007 and 1799, respectively. From Figure 1 of 
Keenan & Norrington (1987), the lower limits imply that n e < 3 x 10 6 cm -3 if T < 60,000 
K. The lines at A4942 and A4989 are contaminated by the strong complex of [O III] AA4959, 
5007, with the result that interesting constraints on T and n e cannot be obtained from these 
lines (Keenan & Norrington 1991). 

Overall, the above results point to a relatively high density for the narrow-line gas, 
which suggests a circumstellar rather than interstellar origin. If the [O III] and [Fe VII] lines 
are assumed to come from gas with similar physical properties, the implied characteristics 
are n e w 2 x 10 6 cm' 3 and T « 20, 000 K, although there is some evidence that the [Fe VII] 
emission comes from a somewhat hotter region. It is likely, because of the low velocity and 
the relatively low pressure, that this gas is the unshocked circumstellar wind. Presumably it 
was photoionized by the UV and X-ray radiation emitted by the supernova right after shock 
breakout. 

One can derive a lower limit to the radial distance of this gas by taking the product of 
the highest observed supernova ejecta velocity and the estimated age. As discussed below, 
the highest velocity is ~ 10 4 km s -1 . The wind electron density in front of the shock then 
becomes 

ne = 5 '° X ( lO- 4 M yr-i ) (iTIm^r (lioW) ^ (1) 

where M is the mass-loss rate from the progenitor star, v w is the wind velocity, and it has 
been assumed that H and He are fully ionized. The reference value of M chosen here is 
close to the highest values that have been observed for Galactic stars and to the strong 
radio emitter SN 1979C (e.g., Lundqvist & Fransson 1988). Either the progenitor star of SN 
1995N was losing mass at an extraordinary rate, or the density of the observed gas is higher 
than the average density in the circumstellar wind. 

The gas density is sufficiently high that the recombination and excitation time scales are 
less than the supernova age, and approximate equilibrium is expected. The ionizing source 
is the X-ray emission, which had an unabsorbed X-ray luminosity of ~ 1 x 10 41 erg s _1 in the 
0.1-10 keV band during July 1996 - January 1998 (Fox et al. 2000). The ASCA observations 
give a temperature estimate of 9 keV if the emission is thermal (Fox et al. 2000). Models 
of X-ray nebulae have been presented by Tarter, Tucker, & Salpeter (1969) and Kallman & 
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McCray (1982). The latter models assume a constant density surrounding a central source 
and generally take T = 10 keV for the photoionizing source, but the calculations can be 
used to approximate the current situation, especially because the [Fe VII] lines are likely to 
be formed in a region optically thin to the ionizing continuum. The ionization properties 
depend on the ionization parameter, which can be defined as £ = L/nr 2 , where L is the 
luminosity of the ionizing source. Taking L — 1 x 10 41 erg s -1 , n = 2 x 10 6 cm' 3 , and 
r = 8.6 x 10 16 cm, we find log£ = 0.8, which is very close to the ionization parameter where 
Fe VII is expected, for a T = 10 keV source (Kallman & McCray 1982). A lower temperature 
for the continuum source would increase the estimated £ for a fixed L. The gas temperature 
of this region in the photoionization models is ~ 15, 000 — 20, 000 K, consistent with the 
observations. This comparison provides support for the overall model and, in particular, for 
the estimated radial distance of the emitting region. 

The intensity of the narrow lines gives another constraint on the physical conditions. 
The hydrogen lines are the best to use because there are no abundance uncertainties. While 
no narrow H I lines can be seen in the day 1007 Keck spectrum because of the strong broad 
component, the higher resolution day 1799 VLT spectrum clearly shows Ha, H[3, and H7. 
For H/9, we find a luminosity of 5.3 x 10 36 erg s _1 at 1799 days. If the emitting region is 
assumed to be at radius r with a thickness Ar and a covering factor /, we have 

t± „ 1.6 x 10- ( - L_ _ ) 3 (j^y 2 , (2) 

where the emissivity of H/3 has been taken for T = 20,000 K (Osterbrock 1989). Although 
the density is uncertain, it is clear that the observations imply a small filling factor. 

Another indication of the small filling factor of the n e = 2 x 10 6 cm~ 3 gas is the 
detection of a substantial X-ray flux by ROSAT. If the gas were pervasive at r = 8.6 x 10 16 
cm, the resulting hydrogen column density, N H « 10 23 cm -2 , would give rise to the complete 
absorption of the soft X-ray flux. The gas is not fully ionized and thus can effectively 
absorb the X-ray emission. In their estimate of the X-ray luminosity, Fox et al. (2000) take 
N H = 0.8 x 10 21 cm -2 based on the expected Galactic absorption. 

Most of the other ions that are observed in SN 1995N require an ionization parameter 
that is similar to that of Fe VII. The ions Fe VI through Fe XI, which are observed, occur over 
a relatively small range of ionization parameter (e.g., Kallman & McCray 1982). However, 
the presence of narrow [O I] emission implies the existence of a separate, lower-ionization 
component. The region could be of low ionization because of an especially high density, but 
it is more likely that the region is neutral because the ionizing photons have been used up 
in the H II region. 
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3.2. Intermediate Component 

The intermediate-width component is best exemplified by the [O III] AA4959, 5007 
(Figure 5), [O II] AA7320, 7330, and [O I] AA6300, 6363 emission lines, which all show box- 
shaped profiles extending to ~ 5, 000 km s _1 . The strongest feature in this component at all 
epochs is the [O III] AA4959, 5007 blend. Other likely members of this component are the 
high-ionization far-UV lines. Mg I] A4571 is a more uncertain member, showing a red deficit 
if it is correctly identified; alternatively, it could be a blend of lines, in particular [Fe II] 
AA4556, 4583. The strongest argument for including Mg I] in this component is its temporal 
evolution, which is similar to that of the [O I] and [O III] lines (§5). There is no indication 
of hydrogen or helium with this type of profile. 

Table 3 lists the fluxes of the lines that we designate as having line profiles coming from 
the intermediate component. As with the narrow lines, the intensity ratios of some of these 
lines provide constraints on the electron density and temperature. 

O III. The [O III] A4959 and A5007 lines are strongly blended with each other, and the 
A4959 line with H/3. A feature at -1, 800 km s" 1 appears in both the A5007 and A4959 lines, 
and the ratio of line strengths appears to be consistent with the 3:1 ratio expected from their 
spontaneous emission probabilities. The red edge of the A5007 line extends to 5, 000 km s _1 , 
but not further. The A5007 line also has a strong feature to the red at +3,000 km s _1 ; 
however, it does not appear in [O III] A4363 or in other lines of this kinematic component, 
so there is a question of whether the feature originates from [O III] or is a blend. Its rapid 
drop in the late-time spectra (see § 5) argues for the blend interpretation. Likely candidates 
for this feature are He I A5049 or [Fe II]. 

Broad [O III] A4363 emission is present, although the blue side of the line is dis- 
turbed by ITy emission. Neglecting the red feature in the A5007 line, we find /(A5007 + 
A4959)//(A4363) « 8.2 from the day 1007 Keck spectrum. In the far-UV, the O III] 
AA1660, 1666 intercombination lines are among the strongest UV lines. Based on the wave- 
length, it is unlikely that the He II A 1640 line contributes much to this feature. We find a 
reddening-corrected ratio [O III] /(AA4959, 5007)//(AA1660, 1666) « 4.0. 

O II. The [O II] AA7320, 7330 doublet is less affected by blending and yields a clearer 
line profile. The feature at —1, 800 km s" 1 is present, but that at +3, 000 km s -1 is not. The 
line profile is asymmetric, although the maximum velocity is 4, 500 km s _1 to both the blue 
and the red. As we show in § 4.1, the asymmetry is mainly caused by blending with He I 
A7281. 

O I. The [O I] AA6300, 6364 doublet is severely affected by the blue wing of Ha. Narrow 
components of the [O I] doublet, as well as [Fe X] A6374.5, also complicate the profile. In § 4 
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we show that the line profile is compatible with that of the other oxygen lines. The auroral 
[O I] A5577 line is likely to be present, although it may be distorted by the corresponding 
night-sky line. We therefore only consider its flux as an upper limit. 

The presence of the O I recombination line at A7774 is of special interest because it 
should be accompanied by strong H I recombination emission with the same line profile, if 
the emission is from material with "cosmic" abundances, since hydrogen is more abundant 
than oxygen by a factor of 10 3 in such material. The absence of the corresponding hydrogen 
emission suggests that the emitting gas is oxygen rich. The O I recombination line at A8446, 
which is expected to have a strength comparable to the A7774 line, is blended with Fe II 
lines, and is also affected by fluorescence (§ 4.2). The emission may be from the oxygen-rich 
ejecta in the supernova and may extend to relatively high velocity because of hydrodynamic 
instabilities with the overlying layer (Fryxell et al. 1991) and/or a high expansion velocity 
of the oxygen core (§ 6.1). 

The lines observed in the short-wavelength G160L HST spectrum are noisy, so there is 
uncertainty in the line identifications. The line widths, ±(3500 — 4000) km s _1 , point to 
memberships of most of the UV emissions with the intermediate kinematic component. The 
line profiles of the strongest lines (C IV, O III, and Si III/ C III) are all consistent with flat- 
topped profiles of width ±3800 km s" 1 . Except for the A1900 feature, which has a central 
peak, the lines appear symmetric within the noise level. The C IV and Si IV lines are also 
affected by interstellar absorption in the Galaxy and the host galaxy. 

By far the strongest line is Lya, whose profile is distorted by geocoronal Lya, and 
possibly by N V A1240. The C IV AA1548, 1551 doublet is clearly seen with a boxy profile. 
The same is true for the feature at A1900, which can be identified as a blend of Si III] A1892 
and C III] AA1907, 1909. The N III] AA1744-1754 lines are clearly present, although their 
flux is uncertain. The N IV] A1486 line is near the limit of detection. 

In Figure 6 we show the density and temperature constraints from all the oxygen- 
line ratios of the HST and Keck day 1007 spectra together. We have used a reddening 
correction E B _ V = 0.11 mag to obtain I([0 I] A7774)/J([0 II] AA7320, 7330) = 0.43, [O III] 
J(AA4959,5007)//(A4363) = 8.2, and [O III] /(AA4959, 5007)//(A1664) = 4.0. The [O I] 
7(A5577)//(AA6300, 6364) ratio is only an upper limit, <0.1, and the curve shown in Figure 
6 therefore only marks an upper boundary to the temperature. Implicit in this analysis 
is that the O I and O III zones are of similar extent, temperature, and electron density. 
While the lines coming from the O II and O III zones are all compatible with T 7000 
K and n e ~ 3 x 10 8 cm~ 3 , the [O I] lines have to arise in a lower temperature and/or 
lower electron density zone. As shown in Chevalier & Fransson (1994, hereinafter CF94), 
both these alternatives are quite natural, because the [O I] lines are expected to arise in a 
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partially neutral zone of lower temperature (see, e.g., CF94 Figures 4 and 5). 

The velocities observed in the intermediate component suggest that this emission is 
from freely expanding ejecta. The ejecta are photoionized by X-rays from the circumstellar 
interaction. This process has been modeled by CF94, and there is general agreement between 
the lines that are predicted to be brightest and those that are observed here. 

3.3. Broad Component 

Because of the strong Ha line, the broad component dominates the optical emission 
from the supernova. The total luminosity in Ha was ~ 2.3 x 10 40 erg s _1 on day 1007, 
using E B -v = 0.11 mag. The line profile of Ha is shown in Figure 5. There is no obvious 
distinction of different components in Ha, as we demonstrate in § 4.1, and no clear evidence 
for any velocity structure within the line, although a slight emission deficit can be seen on 
the red side of the line. The FWHM of the line is ~ 1600 km s _1 on day 1007. The wings 
are very extended and reach at least 10,000 km s _1 . 

Table 4 gives a list of the lines that have profiles similar to that of Ha, although we 
can trace the high-velocity wings above ~ 3000 km s" 1 only for Ha. Besides Ha, the best- 
defined line profile from this component is He I A5876, which is very similar to Ha up to 
~ 3000 km s _1 , where the line is lost in the noise. The broad component is dominated by 
lines of H I, He I, Mg II, and Fe II. 

The interpretation of the Mg II AA2796, 2803 doublet is complicated by absorption 
features arising in gas in our Galaxy and in the host galaxy of SN 1995N. The line peak 
is at negative velocities, at ~ 2791 A, which may be a result of resonance-line scattering 
of the A2796 photons by the A2803 component, as well as the influence of the interstellar 
absorption lines. The line appears considerably broader than Ha, and shows a clear deficit 
of emission to the red of line center (Fig. 5). As discussed in § 4.2, the extended profile is 
likely to be caused by blending with Mg I A2852 and with the strong Fe II resonance lines 
in this region of the spectrum, expected to be present on either side of the wings (multiplets 
399, 391, 380 at ~ 2850 A, and multiplet 373 at ~ 2770 A), due to the Lya fluorescence. 

Although almost all of the many Fe II lines present are blended with each other, it 
appears that they belong to the broad kinematic component. The only reasonably isolated 
Fe II line is A7155, whose profile is well fitted by a scaled Ha profile. In the optical spectra, 
the region 4000-5500 A contains a large number of Fe II lines, often seen in supernovae 
and active galactic nuclei (AGNs; e.g., Filippenko 1989; Osterbrock 1989). Comparing with 
spectra of narrow-line Seyfert 1 galaxies such as I Zw 1 (Phillips 1976, 1977), and the 



- 12 - 



symbiotic star RR Tel (Crawford et al. 1999), we have identified members of Fe II multiplets 
14, 27, 28, 27, 42, 48, 49, and 74 (Table 4). The lines listed in Table 4 are those that can 
be identified with reasonable certainty. Most of the other members of the multiplets are 
blended with strong lines. One should also note that the fluxes of all Fe II lines, with the 
possible exception of the A7155 line, are highly uncertain. 

In the UV, the plateau between ~ 2300 A and ~ 2800 A, except for Mg II, is dominated 
by Fe II resonance lines. Because of the extreme blending of lines, which creates a pseudo- 
continuum, we do not attempt any complete identification in this spectral region. We simply 
note that there are some distinct features such as Fe II AA2614 — 2631 and Fe II AA2506, 2508. 
The latter lines are of special interest, because they arise from cascades from highly excited 
levels, and are usually considered an indicator of fluorescence. Except for this, the UV part 
of the spectrum is similar to that of I Zw 1 (e.g., Fig. 2d in Laor et al. 1997). 

In the near-IR, it is especially interesting to note that many of the features can only 
be identified with Fe II emission based on models of Lya fluorescence of the transitions 
(Johansson & Jordan 1984; Penston 1987; Sigut & Pradhan 1998a). Sigut & Pradhan 
(1998b), who discuss this process in detail for AGNs, provide a detailed list of the predicted 
lines, yielding many identifications for SN 1995N. The best example of Lya fluorescence is 
probably the feature at ~ 9150 A, which is difficult to explain with other mechanisms. Based 
on other Paschen lines, a significant contribution of H I is unlikely, and He I lines are also 
excluded by similar arguments, as well as by their wavelengths. Another candidate, Mg II 
AA9218, 9244, also has the wrong wavelengths. In § 4.2 we show explicitly how the inclusion 
of these lines improves the spectral fit. 

Finally, we point out the very strong He I A10830 line in the VLT day 1769 spectrum, 
which is the only spectrum covering this region. 

Although there are many lines observed with the high- velocity component, the only 
straightforward diagnostic of the physical conditions is given by the lack of forbidden-line 
emission. In the case of the broad-line regions of AGNs, this is taken to indicate a density 
>10 9 cm -3 (e.g., Osterbrock 1989). This should be taken with some caution, however, 
because the gas is only partially ionized by the hard X-rays (CF94); high-ionization lines 
like [O III] AA5007,4959 and C III] A1909 are therefore unlikely to be present. The only 
forbidden lines expected are [O I] AA6300, 6364, which are suppressed relative to permitted 
lines for n e > 10 7 cm~ 3 . As we discuss in § 4.1, the large Ha/H/3 ratio, as well as the Ca II 
line ratios, strongly argue that n e > 10 9 cm" 3 , which is a likely lower limit for the broad 
component. 

The emission is presumably from a radius that is similar to the radius of the narrow-line 
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emission, and the gas is thus exposed to a radiation field that is comparable to that of the 
narrow-line gas. The high density results in an ionization parameter that is a factor >10 3 
smaller than that for the narrow-line gas, or log£ < —2. The result is low ionization of the 
gas, as observed. The ionization is much less than that in the broad-line regions of AGNs, 
where lines like C IV A 1549 are observed. 



To determine abundances as well as deblend the most important lines of the intermediate 
and broad components, we have calculated synthetic spectra based on the combined UV- 
optical HST/Keck spectrum from § 2. The model we employ is a simple two-component 
model, with constant physical parameters and abundances in each of the components. The 
temperature, densities, and abundances are taken as input parameters. For each of the two 
components we assume a given emissivity profile with velocity. The emissivities of the H 
and He lines are given by Case-B recombination, except for Ha, which is multiplied by some 
factor to account for collisions and optical-depth effects. The He I AA3889, 7065 lines are 
also adjusted by a common factor to account for optical-depth effects (see Osterbrock 1989, 
p. 105). All lines from metals are calculated as multilevel atoms for a given temperature and 
density. The ionization fractions for a given ionization stage (such as C III, N III, O III, and 
Si III) are assumed to be the same. The continuum is due to H I, He I, and He II bound-free 
and two-photon continua, consistent with the recombination lines. We do not attempt to 
include the many Fe II lines present, nor lines from the narrow component. The temperature 
and ionization are not calculated self-consistently, as in CF94, so only the relative line fluxes 
for a given ion, or group of ions, are testing the model. The model should therefore mainly 
be viewed as a diagnostic tool. 

For the line profiles we adjust the variation of the emissivity with velocity so that a best 
fit is obtained (eq. 3). For spherical symmetry and for a velocity V(r) = Knax( r /-Rmax), the 
flux at wavelength A is given by 



where S(r) is the source function and r(r) is the optical depth. In the optically thin limit 
S(r)[l — exp(— t)] ~ j(r), where j is the emissivity. The limits of integration are given by 
x = ct(X — A )/A , and R = V max t, where V max is the maximum ejecta velocity and c is the 
speed of light. This characterization of the line profile is appropriate for freely expanding 
ejecta; free expansion is assumed in our synthetic spectral models. As discussed in § 6.1.3, 
this assumption may not apply to the broad component. In this case, the line profile should 
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be viewed as an empirical fit that is applied to the set of lines. 

4.1. Results 

We find that the profile of the H I, He I, Mg II, and Fe II lines can be well reproduced 
with an emissivity [or, if r 3> 1, a source function S(r)] which is constant inside 1000 km s" 1 
and decreasing like j(V) oc V' 4 - 6 outside, out to about 15,000 km s -1 . For the metal lines, 
except for Mg II and Fe II, we have used a const ant-emissivity shell with j = constant for 
2500 < V < 5000 km s _1 . 

In Figure 7 we show the spectrum of a model where all elements, including H, He, and 
metals, have the same density and temperature, 4 x 10 6 cm~ 3 and 14, 000 K. The total 
metallicity is close to solar. Note that the temperature, density, and metallicity are only 
indicative, and depend on the assumed total emissivity, f n e nr 2 dr. The fit gives a good 
general representation of most lines, in terms of both the relative fluxes and line profiles. In 
particular, many of the profiles that seem fairly complicated at first sight are well reproduced 
as a result of blends by several lines. Good examples are the [O I] + Ha complex, and the 
[O III] + H/9 feature at 4800-5000 A. The model, however, severely underproduces the 
strength of the important O I A7774 recombination line. 

This could have been anticipated already from Figure 6, which shows that the relative [O 
III] ratios can be reproduced either with a high-temperature, low-density model, as in Figure 
7, or with a low temperature and a high density. The [O II] AA7320 + 7330/O I A7774 ratio 
can, however, only be reproduced with the latter temperature/density combination. The 
relative fluxes of the oxygen and hydrogen/helium lines then require a very high oxygen 
abundance, unless the O II and O III zones are unreasonably large compared to the zone 
responsible for the H I and He I recombination lines. For the same extent, the relative flux 
in O I A7774 and H I lines is given by [j cff (0 I A7774)n(0 II)]/[j cff (H I)n(H II)]. Because 
jes(0 I A 7774) ^ 2 x j cS (H.{3), this shows that n(0 II) and n(H II) must be of the same 
order of magnitude. 

As an alternative, we show a model with a greatly increased metal abundance (Fig. 8). 
Although these lines are likely to arise in different ionization zones, we find that we can get 
an acceptable fit to all lines with a density ~ 3 x 10 s cm~ 3 and a temperature ~ 7, 000 
K. All metal lines, except Mg II and Fe II, are assumed to have a shell-like emissivity 
profile. Compared to Figure 7 we see that, because of the high oxygen abundance, now both 
[O II] AA7320, 7330 and O I A7774, as well as the strengths relative to the H I lines, are 
well reproduced. In addition, the O I A8446 recombination line contributes appreciably to 
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the feature at ~ 8450 A. We also note the good fit to the O I A7774 line profile, further 
strengthening this identification. 

To reproduce the observed Ha/H/5 ratio, the Ha strength has to be increased by a factor 
of ~ 8 above recombination. The reddening-corrected Lya/Ha ratio is ~ 2.3, while the Case- 
B ratio is ~ 24. The flux of Lya is, however, highly uncertain, because the line profile is 
disturbed by geocoronal Lya, in spite of the recession velocity of SN 1995N. As the self- 
consistent models in CF94 show, this occurs because of collisional excitation in combination 
with a high optical depth for Ha in the partially ionized zone. As both the calculations 
in CF94 and similar models for AGNs (e.g., Kwan & Krolik 1981) show, this requires the 
presence of hard X-rays from the shocked region in combination with a density >10 8 cm -3 . 
The higher-order Balmer lines are well explained by Case-B recombination. 

The helium lines are well reproduced by pure recombination, except for He I A3889 and 
He I A7065. The former is down by a factor of ~ 0.5, and consequently the He I A7065 line is 
increased by a factor of 2.0 x 7065/3889 3.6, which strongly indicates that the He I A3889 
line is optically thick (e.g., Osterbrock 1989, p. 105). 

The relative ratio of the He I and H I fluxes results in a He II/H II abundance of ~ 0.35 
by number. Because the He II and H II ionization zones are likely to coincide (e.g., CF94), 
the He II/H II ratio is probably close to the elemental He/H abundance ratio. Strong He I 
lines are also present in SN 1986J (Rupen et al. 1987; Leibundgut et al. 1991) and at early 
times in many SNe Iln (e.g., Filippenko 1997). 

The jump at A ~ 3700 A is mainly caused by high-order Balmer lines and the Balmer 
continuum. The two-photon continuum dominates in most of the UV, as well as in the 
optical. The form of the Balmer continuum is sensitive to the temperature in the zone where 
this is formed, and a somewhat better fit is obtained for T « 9000 K. In a one-zone model, 
however, this temperature gives a worse fit to the metal lines. More detailed models (CF94) 
show that both temperature and ionization vary within the emitting regions. 

The Ca II lines in the spectrum are interesting. The IR triplet lines at AA8498, 8542, 
8662 are clearly present, while for the H&K lines and the forbidden AA7291.5, 7323.9 lines we 
can only set upper limits. Fransson & Chevalier (1989) and Ferland & Persson (1989) have 
discussed the relative strengths of these lines as a function of density and temperature. Using 
their results we find that the fact that /(AA7291.5, 7323.9)//(AA8498, 8542, 8662) < 0.25 
and J(H&K)//(AA8498, 8542, 8662) < 0.12 means that the density must be >10 9 cm" 3 . 
This result is only weakly sensitive to the temperature. 

The parts of the spectrum which are not well fitted by the model can in most cases be 
identified with regions of strong Fe II emission. This includes the UV between 2300 — 3000 
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A, some regions in the range 4200 A to 5300 A, the A7155 line, and the fluorescence lines 
above ~ 8000 A (see below). 

In Figure 9 we show the line profiles and fits of the Ha, [O I], [O II], and [O III] lines 
to illustrate the fits for our assumed radial emissivity model. The figure also demonstrates 
the complex spectrum which results from the blending of the individually fairly simple line 
profiles. 

Turning now to the other metal lines coming from the intermediate- velocity component, 
the UV line-intensity ratios require N/C = 2.5 — 5 and O/N = 3.5 — 4.7, by number. The high 
N/C ratio, in combination with the relatively large O/N ratio, strongly indicates incomplete 
CNO processing in this component. 

While the profiles of the high-ionization C, N, O, and Si lines in the far-UV are con- 
sistent with the same shell-like emissivity as the oxygen lines, this emissivity model shows 
a bad fit of the Mg II AA2796, 2803 line. The peaked profile shows that this line is instead 
likely to originate in the H and He-rich broad-line component, rather than in the metal-rich 
component. 

4.2. Fe II — Lya fluorescence 

Based on the Fe II - Lya fluorescence calculation by Sigut & Pradhan (1998a,b), we have 
modeled the near-IR spectrum of SN 1995N. We have used the same model as in Figure 8, 
but now added the Fe II lines, with wavelengths and fluxes from the list of Sigut & Pradhan. 
The line profiles were assumed to be the same as that of Ha. In the upper panel of Figure 
10 we show the resulting spectrum, while in the lower panel we show for comparison the 
spectrum without Fe II lines. The observations are from the day 1769 VLT spectrum, but 
the day 716 and 1037 Lick spectra show the same features. The day 1007 Keck spectrum 
only extends to ~ 9000 A. 

As can be seen, there is a striking improvement with the observed features between 
8000 A and 9500 A when the Fe II lines are included. In particular, the feature at ~ 9150 A 
can be identified with the Fe II AA9071, 9128, 9177 lines, and provides a good fit both in 
wavelength and profile of the blend. 

One of the strongest expected fluorescence-pumped Fe II lines is at A8451, and there is 
indeed a strong feature at this wavelength. Another possible identification for this feature is 
the recombination line of O I A8446. This line is normally accompanied by O I A7774, which 
is not present in the broad component (although it is visible in the intermediate component; 
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see § 4.1). The recombination contribution to the A8446 line is expected to be close to the 
A7774 line, and is therefore insufficient to explain this strong feature. The A8446 line can 
be produced by fluorescence through the Ly/3 line (Grandi 1980). Leibundgut et al. (1991) 
identified the O 1 A8446 line in the spectrum of SN 1986J. The Ly/3 -01 A8446 fluorescence 
also results in a strengthened O I A1300 multiplet. Due to its narrowness, however, this is 
likely to be absorbed by the O I in the host galaxy. In our Keck day 1007 spectrum of SN 
1995N, the center of the emission feature gives a better match to Fe II A8451 than to O I 
A8446. An additional argument in favor of this mechanism is that Fe II AA2506, 2508, which 
is indicative of Lya pumping, is the strongest Fe II feature in the UV. This has previously 
been seen in rj Carina, RR Tel, and KQ Pup (Baratta, Cassatella, & Viotti 1995; Johansson 
& Jordan 1984; Redfors & Johansson 2000). Other lines in the UV expected to be enhanced 
by fluorescence are the Fe II multiplets 399, 391, and 380 at 2850 A, multiplet 373 at 2770 A, 
and multiplet 363 at 2530 A. These most likely contribute to the broad wings of the Mg II 
lines, but are too blended to be uniquely identified. IR lines at ~ 1.7 fim are also expected 
to be strong (Johansson & Hamann 1993), but are outside our observed wavelength range. 

Given that the model was intended for conditions in an AGN, the fit is surprisingly 
good. Only at >9300 A are there a few features (such as at ~ 9573 A) which are not well 
reproduced, but the noise level at these wavelengths is considerable. One should also note 
that the relative strengths depend on the uncertain atomic data, as well as on the exact Lya 
profile. In particular, the larger line width in SN 1995N compared with those assumed in the 
AGN calculation by Sigut & Pradhan (1998a,b) model can change the results appreciably. 
Sigut & Pradhan assumed a static Lya profile with only damping included. An expanding 
medium means that the Lya line can feed many more fluorescence channels, which is likely 
to enhance the fluorescence effects. 

Finally, we note the strong He I A10830 line in the day 1769 spectrum. Although the 
region beyond ~ 9300 A becomes increasingly noisy, this feature clearly stands out above 
the noise level. Compared to the strength expected by pure recombination, it is enhanced by 
collisional excitation by a factor of ~ 8, which is common in many objects (e.g., Osterbrock 
1989, p. 111). 

5. TIME EVOLUTION 

In Figure 11 we show the evolution of the reddening-corrected Ha luminosity with time 
on a linear-log scale. We see that the evolution can be well approximated with an exponential 
decay, L(Ha) oc exp(— t/r), where r = 460 days. This is much longer than the 56 Co decay 
time scale, and we conclude that the circumstellar interaction is responsible for the longevity, 
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although there is no obvious explanation for this particular time dependence. 

An efficient way of displaying the change of the spectrum with time is to plot the ratio 
of two spectra as a function of wavelength. In Figure 12 we show a running average taken 
over 2 A of the day 1799 VLT and day 1007 Keck spectra. The presence of two separate 
components in the spectrum is apparent. While most of the spectrum has decreased by a 
factor of ~ 6, the intermediate-component [O I], [O II], and [O III] lines only decreased by a 
factor of 2-3. The O I A7774 line shows a similar decrease, while the A8450 feature follows 
the H I lines, consistent with it being mainly due to Fe II fluorescence lines. In addition, 
Mg I] A4571 belongs to the class of slowly declining lines. 

Relative to H/5, the peak fluxes of the [O I] and [O II] lines increased by a factor of 
~ 3 between 1997 and 1999, while the [O III] line increased by a more modest factor of 
~ 1.6 between these observations. The [O III] AA4959, 5007/A4363 ratio remained nearly 
constant at ~ 7.9, while the O I A7774/[0 II] AA7320, 7330 ratio decreased from ~ 0.43 
to ~ 0.20. Repeating the analysis in § 3.2 (Fig. 6), we find that the temperature in 
the oxygen-rich gas had not changed much, while the electron density decreased from ~ 
3 x 10 s cm" 3 to ~ 1 x 10 s cm" 3 . This is less than expected from a simple n e oc t~ 3 law, 
but our analysis assumes that the relative sizes of the O II and O III zones remain constant, 
and that the temperatures and densities in these two zones are similar. Neither of these 
conditions is obviously satisfied. In addition, the supernova may be older than our estimate, 
and the uncertainties in the relative fluxes are considerable. Nevertheless, the behavior is 
qualitatively that expected for expanding ejecta. 

In Figure 13 we show the relative change of Ha, [O I] AA6300, 6364, [O II] AA7320, 
7330, and [O III] AA4959, 5007 on a velocity scale. While the [O I] and [O II] line profiles 
changed by roughly the same factor over the whole velocity range, —5000 to +5000 km s _1 , 
the red side of the [O III] line decreased by a large factor relative to the blue. The most 
conceivable reason for the rapid [O III] change is that blending with another line, decreasing 
similarly to the H I and He I component, causes the apparent decrease in the red wing. The 
most likely candidates for this are the He I AA5015, 5049 and Fe II lines. The blue side 
behaved in a manner similar to the [O I] and [O II] lines. We therefore conclude that there 
is little change in the overall profiles of the oxygen lines. When comparing the width of the 
best-isolated [O II] line, we find a nearly constant blue extension, while the width on the red 
side decreased by ~ 450 km s _1 during this period. 

The evolution of the Ha line profile differs dramatically from that of the oxygen lines. 
As shown in the upper-left panel of Figure 13, the red wing up to ~ 3000 km s _1 decreased 
by almost a factor of 4 relative to the blue, with a maximum decrease at ~ 1100 km s _1 . 
To further illustrate the different evolution of the wings of Ha, we have reflected the blue 



-19- 



component in the spectra from 1996 - 1999 (Fig. 14), so that the blue and red wings can be 
compared directly. The fluxes have been normalized to the same peak value. We see that the 
day 716 Lick spectrum (dashed lines) showed nearly identical red and blue wings. The day 
1007 Keck spectrum (solid lines) shows a clear evolution in that for a given velocity the red 
wing decreased faster than the blue wing. This is even more pronounced in the 1999 VLT 
spectrum, where at ~ 1500 km s _1 the flux of the red wing is only ~ 20% that of the blue 
wing at the same velocity. Equivalently, for a given intensity the red wing has a much smaller 
extension. For example, at ~ 10% of the peak flux the blue wing extends to ~ 1900 km s _1 , 
while the red wing only extends to ~ 1080 km s _1 . Relative to the peak flux, the blue wing 
shows a steady increase in flux with time, while the red wing displays the opposite trend. 
Similar evolution is seen for both H/3 and He I A5876, although the continuum contribution 
and line blending make it difficult to compare directly to Ha. Baird et al. (1998) have 
presented evidence that the line slowly narrows with time. We can therefore confirm this 
trend to the latest epochs. The Ha/H/3 ratio in the 1999 spectra decreased somewhat relative 
to the 1997 and 1996 spectra. This is consistent with the decreasing density during this time 
interval. 

The intensity ratio of the He I lines to the Balmer lines does not change much. The two 
He I lines least affected by blending, He I AA5876, 7065, have almost the same ratio relative 
to B./3 at all epochs, and are all reproduced with a relative He II/H II abundance of ~ 0.35 
by number. Together with the similar line profiles, this argues for the fact that the He I 
and Balmer lines arise in the same region. Their relative flux should therefore reflect their 
abundance ratio. 

In most cases the fluxes of the narrow lines are uncertain, especially in the day 1799 
spectrum. For the strongest lines, like [O III] A5007 and [Fe VII] A6087, the decrease is 
by a factor of 5-10, with no systematic distinction between lines from different elements. 
Although the flux of the [O III] A4363 line is especially uncertain due to the underlying 
broad component, there is an indication for an increase in the [O III] A5007/A4363 ratio 
from ~ 4.1 to ~ 6.6. However, because this ratio depends on both density and temperature, 
one cannot directly translate this into a unique change in either of these parameters. 

6. DISCUSSION 

6.1. Nature of the Emission Regions 

The general scenario for the emission from SN 1995N and related supernovae (generally 
Type Iln) is interaction of the supernova ejecta with circumstellar gas. Our spectra of SN 
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1995N show that the interaction is complex, but that there is sufficient diagnostic information 
on the emission components to deduce some of the interaction details. All the emission 
components are likely to be heated and ionized by the energetic radiation from the interaction 
region. 

6.1.1. Narrow Lines 

Based on the unresolved line widths, the density of ~ 2 x 10 6 cm -3 , and the low filling 
factor, the narrow-line emission can be identified with clumps in the preshock circumstellar 
gas. The velocity is likely to be < 50 km s _1 , typical of the circumstellar media of red 
supergiants. 

6.1.2. Intermediate Component 

The intermediate-velocity component is identified with unshocked supernova ejecta. 
The maximum velocity in this component suggests either that the reverse shock front has 
a minimum velocity of ~ 5,000 km s -1 , or alternatively that the emission is coming from 
a component which is photoionized by the shock but only extends in abundance to this 
velocity. The fact that this component is only seen as dense, oxygen-rich gas makes it natural 
to identify it as processed gas from the core. Assuming that the oxygen core expands at 
5, 000 km s^ 1 , the typical oxygen density will be 

n ~ 4.1 x lO 5 /- 1 ( ( V T e ; ) 3 ( ^— 1 3 cm" 3 , (4) 

J \2 Mq) V5000 km s" 1 J V 1000 da Y s / 

where / is the filling factor of this component in the core, and we have assumed an oxygen- 
dominated composition. With / ps 0.1, as was the case for SN 1987A (e.g., Spyromilio & 
Pinto 1991), the density is of the same order as found from the [O III] and similar lines. 
We note that the density deduced from the models is based on the assumption that the 
temperatures and densities are the same in the O I, O II, and O III zones. 

A velocity of ~ 5, 000 km s _1 is high for the oxygen core of a Type IIP or Type IIL 
supernova, but in the range found in Type Ib/Ic supernovae, as well as in the Type lib 
SN 1993 J (e.g., Houck & Fransson 1996; Matheson et al. 2000). This is an important result 
because it is consistent with a picture where the progenitor has lost most of its hydrogen-rich 
envelope before the explosion. 
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6.1.3. Broad Component 

While most of the hydrogen-rich gas comes from gas moving at velocities of <2000 km s _1 
there is a smooth extension to ~ 10 4 km s" 1 . Any model must explain the wide range of 
velocities for this component. In particular, the fact that there appears to be hydrogen-rich, 
cool gas moving at up to 10, 000 km s" 1 implies that there are some regions where the shock 
region extends considerably beyond the reverse shock front. 

To explain this, one may identify three main alternatives for the broad, hydrogen-rich 
component. The first possibility is that it comes primarily from unshocked, freely expanding 
ejecta, extending from low-velocity hydrogen in the core up to the maximum ejecta velocity. 
A second possibility is that it comes from shocked ejecta gas, behind the reverse shock. The 
final possibility is that it comes from shocked circumstellar gas, either in the form of clumps 
or a more uniform medium. Our synthetic spectral models in § 4 are appropriate for the 
first possibility. 

Before discussing these possibilities, we note that the Ca II lines, the absence of forbidden 
lines, and the Ha/H/3 ratio indicate a density >10 9 cm" 3 for this component (§ 4.1). A 
further independent lower limit to the density can be obtained by considering the state of 
ionization in this gas. The luminosity in Ha was ~ 2.3 x 10 40 erg s -1 on day 1007 (§ 3.3). 
Typically, the conversion efficiency from X-rays to Ha is ~ 1% (CF94), indicating a total 
X-ray luminosity of 2 x 10 42 erg s _1 . This is much more than was observed by ROSAT (Fox 
et al. 2000), but as calculations for SN 1993 J showed (Fransson, Lundqvist, & Chevalier 
1996), almost all the X-ray emission can be absorbed in the ejecta and shocked gas, and only 
the hard X-rays escape the supernova. With this luminosity the ionization parameter on day 
1007 is £ «3 1 x 10 2 (n/10 6 cm" 3 )" 1 . For a uniform density, hydrogen-dominated envelope 
expanding at ~ 10, 000 km s _1 , the density is 

a, ^ f -i( M ™\( V ™ VV 1 V 3 -3 ^ 
n w 4.1 x 10 f cm . (5) 

J \i m q J Vio,ooo km s - v V 1000 da y s / 

Therefore, if the emission came from unshocked ejecta and / > 0.1, the ionization parameter 
would be very high, £ > 10, and one would expect strong emission from ions like [O III] 
from this component, contrary to the observations. Instead, it is dominated by Mg II and 
Fe II emission, in addition to the H I and He I lines, characteristic of a partially ionized, 
high-density region (e.g., CF94). It is therefore unlikely to come from a uniform-ejecta 
component, and only if clumping of this component is very high (/ < 0.01) is an ejecta 
origin likely. We therefore think the first possibility given above is unlikely. 

The density in the broad component is much higher than that in the narrow-line com- 
ponent, but the temperature is expected to be only slightly smaller. The conclusion is that 
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the pressure in the gas is relatively high and one possibility is shocked circumstellar or ejecta 
gas. The gas pressure is p/k = 2 x 10 13 ngT 4 cm -3 K, where ng is the H density in units of 
10 9 cm" 3 and T 4 is the temperature in units of 10 4 K. The pressure generated by a shock 
front moving into gas with a density of 10 6 cm" 3 is p/k = p Q v 2 s /k = 1.5 x 10 14 w^ 3 cm" 3 K, 
where v s s is the shock velocity in units of 10 3 km s _1 . If the emission is from shocks moving 
into a medium with a density of this order, the required high density can easily be produced. 

We can identify two possibilities for such a shocked gas: it could be either clumps caused 
by instabilities in the shocked ejecta behind the reverse shock, or dumpiness and instabilities 
of the shocked circumstellar medium. This is similar to the two scenarios discussed by Chugai 
& Danziger (1994) for SN 1988Z. 

Shocked ejecta: This is a version of the scenario discussed for "normal" interacting 
supernovae like SN 1979C or SN 1993 J (CF94; Fransson et al. 1996). The fact that the line 
profiles in SN 1995N are peaked as opposed to the boxy profiles seen in SN 1993J, however, 
makes it necessary to modify this scenario. One possibility is that the circumstellar medium 
is highly asymmetric, as in the picture discussed by Blondin, Lundqvist, & Chevalier (1996). 
Indications for a highly asymmetric circumstellar medium are seen in some red supergiants 
in a superwind phase (see § 6.2). Most of the hydrogen emission could in this case come from 
the cool, dense shell expected behind the reverse shock (CF94). The shape of the Ha line 
profile then reflects the decreasing emission measure as one moves away from the equator. 
The maximum ejecta velocity is likely to occur at the poles, while the strongest interaction, 
with most of the energy release, occurs in the high-density circumstellar gas at the equator. 
The deceleration of the ejecta will therefore be most rapid here. If the hydrogen envelope 
mass is low, the reverse shock at the equator may rapidly recede to the massive oxygen 
core. This could qualitatively explain the presence of the high-ionization, oxygen-rich shell 
observed in the intermediate-component line profiles. The hydrogen-rich gas in the cool shell 
behind the reverse shock could then be excited by the X-rays from the shock, giving rise to 
an extended partially ionized zone. A prediction is that the maximum velocity of the shocks, 
and therefore the line widths, should decrease with time. 

This general scenario for SN 1995N receives support from observations of the related 
supernova SN 1986 J in NGC 891, which had a strong Ha line with FWHM « 700 km s" 1 
and oxygen lines with widths of several 10 3 km s" 1 (Leibundgut et al. 1991). It was a 
strong radio source, and VLBI images showed a shell with protrusions extending at least a 
factor of two beyond the shell radius (Bartel et al. 1989, 1991). Images on different dates 
revealed expansion; the outer protrusions are moving at ~ 15, 000 km s" 1 , while the shell is 
expanding at ~ 5,000 km s" 1 . These images suggest a picture that is consistent with that 
deduced from the spectra of SN 1995N. 
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The origin of the protrusions is uncertain, but there are computational simulations 
that suggest that they could arise from asymmetries or inhomogeneities in the circumstellar 
medium. Blondin et al. (1996) found that supernova expansion into a circumstellar medium 
with a density deficit in the polar direction could produce polar protrusions with a radius 2-4 
times that of the main shell. The problem in comparison to the images of SN 1986 J is that 
there are more than two protrusions. Jun, Jones, & Norman (1996) carried out a simulation 
of the interaction of a supernova with a clumpy medium and found that the Rayleigh- Taylor 
instability that develops where the ejecta are decelerated by the surrounding medium is 
accentuated by the presence of clumps. The Rayleigh- Taylor fingers reach the outer shock 
front. The protrusions in the simulation did not reach the extent required for SN 1986J, 
but the simulation was run for a limited amount of time and the density of the surrounding 
medium was uniform. The circumstellar wind case remains to be investigated. 

Shocked circumstellar clouds: In this scenario, the emission comes from the cool gas 
behind the slow radiative shocks propagating into the same clumpy gas clouds as are re- 
sponsible for the narrow-line emission. The blast wave may then propagate in the intercloud 
medium, while the high-density circumstellar clumps are left behind. The maximum ejecta 
velocity will decrease only slowly, while the shocked clouds will be accelerated on a time 
scale which depends on the column density of the clouds. As simulations show (e.g., Klein, 
McKee, & Colella 1994; Jun et al. 1996), however, the cloud may gradually be evaporated, 
and finally be dissolved and merge with the ejecta gas. In this scenario one expects a smooth 
distribution of gas, with velocities ranging from newly shocked clouds with low velocity, or 
clouds with very high column density, to gas which has been accelerated to the ejecta ve- 
locity. This agrees qualitatively with the observed Ha line profile. The observed gradual 
increase of the flux in the blue wing with time is also consistent with this scenario. The 
decrease in the red wing, however, is not naturally explained, unless it is caused by dust 
absorption, either from pre-existing dust in the circumstellar gas, or from newly formed dust 
in the ejecta, as discussed for SN 1998S by Gerardy et al. (2000). 



6.1.4- Location of the Energy Source 

An important question is the location of the energy source for the excitation of the 
gas, and there are several possibilities. First, the source could be the shock separating the 
high-velocity hydrogen ejecta from the circumstellar gas. In the uniform-medium case this 
will be at the highest velocity and also maximum ejecta radius. This is the model discussed 
by CF94, and applied to SN 1993J by Fransson et al. (1996). However, it is difficult to 
understand the presence of the high-ionization, oxygen-rich shell at only ~ 5000 km s -1 , 
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since most of the X-rays will be absorbed by the unprocessed ejecta between the shock and 
the oxygen core. 

In the clumped circumstellar medium case the bow shocks may be distributed over 
a large radius, and immersed in the expanding ejecta, depending on the time of shocking 
and their column densities. The oxygen-shell emission is again difficult to explain. It may, 
however, explain the distributed nature of the hydrogen emission. 

6.2. Type Iln Progenitors 

Some clues to realistic scenarios for the structure of the circumstellar medium can be 
obtained from observations of Galactic supergiants in very advanced evolutionary stages. 
While most red supergiants, like a Ori, have fairly modest winds with mass-loss rates 
<1(T 5 M yr-\ there is a handful of massive stars having; more extreme mass-loss rates, 
10~ 4 — 10~ 3 M yr _1 , often referred to as "superwinds." Van Loon et al. (1999) find that 
two out of eight red supergiants in an LMC sample have mass-loss rates ^,10 -3 M Q yr _1 , 
while the others have moderate mass-loss rates ^,10~ 5 M yr _1 . 

The best-studied superwind cases are VY CMa and IRC +10420. The nature of the 
superwind phase is not yet clear, but can be related either to a specific evolutionary phase, 
occurring for most red supergiants, or to a phase only occurring for a specific mass range. 
A hint that the latter is the case comes from the fact that both VY CMa and IRC +10420 
are very luminous stars, with L (2 — 5) x 10 5 L . From evolutionary models Wittkowski, 
Langer, & Weigelt (1998) estimate a zero-age main sequence (ZAMS) mass of 30-40 M Q 
for VY CMa. The mass-loss rate is determined to be ~ (2 — 3) x 10 -4 M yr _1 for a wind 
velocity of 39 km s" 1 (Danchi et al. 1994). 

Several authors have found an asymmetric distribution of the circumstellar gas and dust 
around VY CMa. From HST imaging Kastner & Weintraub (1998) find an asymmetric core 
region with wavelength dependent FWHM « 10 15 cm in the optical. The star itself has a 
radius of ~ 1.2 x 10 14 cm. Outside of this is an asymmetric reflection nebula with an extent 
of ~ 6 x 10 16 cm, with the core close to its edge. Based on symmetry arguments, Kastner 
& Weintraub speculate that this is only one half of a bipolar nebula, with the other half 
obscured by an opaque equatorial dust disk. The same conclusion is reached by Monnier et 
al. (1999) from high-resolution observations of the core and extended structure. 

The supergiant IRC +10420 is similar to VY CMa in terms of mass- loss rate and mass, 
but probably in an even later evolutionary state, intermediate between a red supergiant and 
a Wolf-Rayet star. There is also evidence for a disk-like geometry in this case, based on line 
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profiles (Jones et al. 1993) and on HST and IR imaging (Humphreys et al. 1997). 

Wittkowski et al. (1998) propose that the asymmetric mass loss arises as a result of 
rapid rotation in the phase preceding the evolution off the Hayashi line when the major part 
of the envelope is lost. VY CMa is now estimated to have a mass of ~ 15 M , meaning 
that it has lost ~ 25 M to the circumstellar medium. The helium fraction at this phase is 
~ 0.4 by mass. This is lower than, although qualitatively consistent with, our results for SN 
1995N, but is obviously sensitive to mixing prescriptions, the exact evolutionary status, and 
other factors. As an alternative, Soker (2000) proposes that the asymmetric distribution is 
caused by the concentration of magnetic spots in the equatorial plane, which through the 
lower photospheric temperature causes dust formation preferentially at the equator. 

Based on these cases, it seems probable that SN 1995N, and other Type Iln supernovae, 
arise from progenitors similar to VY CMa or IRC +10420. This is consistent with the high 
mass-loss rates needed to explain the strong interaction, as reflected in the total luminosity. 
A low hydrogen envelope mass is needed to explain the large oxygen core velocity. The 
observed high He/H ratio is also consistent with high mass loss. An asymmetric circumstellar 
distribution will result in faster expansion of the supernova ejecta in the polar directions, 
while it will be slowed down by the equatorial disk or torus as modeled by Blondin et al. 
(1996). The radius of the inner boundary of the dusty disk in VY CMa is estimated to be 
~ 10 15 cm (Monnier et al. 1999). If the density in the disk is large enough, the radius of 
the inner edge would not be much affected by the supernova ejecta. The relative size of the 
ejecta and obscuring disk would therefore increase with time. If the disk is opaque due to 
dust, an increasingly large fraction of the ejecta on the far side of the supernova would be 
obscured. This could explain the fact that the red side of the Ha line decreased considerably 
faster than the blue side. The exact temporal evolution depends on the degree of asymmetry 
of the ejecta. 

Evidence for dusty environments comes from several Type II supernovae. Especially 
interesting are the observations of SN 1998S, where Gerardy et al. (2000) find strong evidence 
for dust emission. As Gerardy et al. discuss, although this dust could in principle be formed 
in the ejecta, perhaps a more likely interpretation is pre-existing circumstellar dust, heated 
by the supernova (see also Fassia et al. 2001). 

From observations of H 2 masers in VY CMa, Richards, Yates, & Cohen (1998) estimate 
a density of ~ 5 x 10 9 cm" 3 in the masing regions, while other species give a density of 
~ 10 8 cm" 3 , indicating a considerable dumpiness of the circumstellar medium. One then 
speculates that these masing regions could be related to the high-density clumps needed to 
decelerate the supernova blast wave. 
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Finally, the small fraction of supernovae in this class, <5% (Cappellaro et al. 1993), is at 
least qualitatively consistent with them coming from a fairly small fraction of the high-mass 
stars. 



7. CONCLUSIONS 

Late-time spectral observations of SN 1995N are presented, covering the UV to near-IR 
regions. Three kinematic components are found, which most likely are powered by X-rays 
from the interaction of the ejecta and the circumstellar medium of the progenitor. Spectral 
modeling shows evidence of an increased helium abundance, as well as CNO burning. The 
first evidence for Lja pumped fluorescence of Fe II lines in supernovae is found. It is likely 
that other supernovae interacting with their circumstellar media show similar signatures, 
especially in the near-IR. While the narrow lines come from unshocked circumstellar gas and 
the intermediate- velocity component from processed ejecta, the high-velocity unprocessed gas 
has a more uncertain origin. Different possibilities based on either clumpy or asymmetric 
circumstellar media are discussed, without strongly favoring either of these possibilities. 
The temporal evolution of the processed and unprocessed gas is distinctly different, again 
emphasizing the different origins of these components. We propose that the progenitors 
of the narrow-line supernovae (Type Iln) are similar to red supergiants in their superwind 
phase, as most of the hydrogen-rich gas is expelled in the last ~ 10 4 years before explosion. 
This explains the high oxygen core velocity. We also present evidence that the circumstellar 
medium is either asymmetric or clumpy. 
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Fig. 2.— HST PC F814W image of SN 1995N from 22 July 2000. The SN has I « 20.5 
mag. 
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Fig. 3. — Spectra of SN 1995N within 5 years after discovery. All except the 11 April 1997 
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Fig. 4.— Combined day 943 (HST) and day 1007 (Keck) spectra of SN 1995N, with iden- 
tifications of the strongest lines. We have omitted identifications of the many Fe II lines in 
the UV and visible ranges. Wavelengths are in the rest frame of the supernova. 
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Fig. 6. — Temperature vs. density plot for the oxygen lines. The [O I] AA6300, 6364/A5577 
curve marks only an upper limit to the temperature for a given electron density. 
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Fig. 7. — Synthetic spectrum superimposed on the day 1007-1037 spectra, corrected for 
reddening. The model includes H, He, C, N, O, Mg, Si, and Ca, but no Fe II lines. The 
metallicity is close to solar, T = 14000 K, and n e = 4 x 10 6 cm" 3 . Note the severe under- 
production of O I A7774, caused by the normal metallic abundances. 
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Fig. 8. — Same as Fig. 4, but now with T = 7500 K, n e = 3 x 10 s cm -3 , and a high metal 
abundance. The O I A7774 line is now well-reproduced relative to the other O I and O II 
lines. 
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Fig. 9. — Spectral fit and decomposition of the Ha, [O I], [O II], and [O III] line profiles from 
the day 1007 spectrum. The velocity scale is counted from Ha, [O III] A5007, and [O II] 
A7325, respectively. The dashed lines give the most important separate contributions to the 
total line profile. The narrow components of the [O I], [O II], and [O III] lines are apparent 
as the excess emission. 
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Fig. 10. — The observed near-IR region of the spectrum of SN 1995N on day 1769, together 
with the synthetic spectrum expected based on the Lya fluorescence mechanism (thick line). 
The upper panel includes the Lya fluorescence, while the lower does not. Note especially 
the greatly improved fits at ~ 8500 A and ~ 9150 A. Also note the strong He I A10830 line. 



-41 - 



41 



^40.5 



TO 
W 

SB 
CD 



B 40 
o 



39.5 - 



39 



h 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



X • 

s 

s 



J I I I I I I I I I I I I I I I I I L 







500 



1000 
Epoch (days) 



1500 



2000 



Fig. 11. — Evolution of the Ha luminosity with time. The dashed line gives a least-squares 
fit to the data. 
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Fig. 12. — Ratio of the day 1799 and day 1007 spectra with identifications of the most 
important lines. Note the peaks of the [O I], [O II], [O III], and Mg I lines, showing that 
these belong to a component separate from the H I, He I, and Fe II lines. 
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Fig. 14. — Evolution of the Ha line profile from 1996 (Lick; dashed line), to 1997 (Keck; 
solid line), to 1999 (VLT; dotted line). Note the strong decrease in the red wing of the line 
relative to the blue wing in the 1997 and 1999 spectra. 
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Table 1. Log of observations 



UT Date 


Epoch a 


Telescope 


Wavelength Range 


Resolution 




days 




(A) 


(A) 


24 May 1995 


321 


Lick 3 m 


3120 - 10490 


6-11 


6-7 July 1995 


364-365 


Lick 3 m 


3120 - 10490 


6-11 


1 August 1995 


390 


Lick 3 m 


4260 - 7040 


6 


22 June 1996 


716 


Lick 3 m 


3120 - 9880 


6 


6 February 1997 


943 


HST 


1140 - 4780 


2-7 


11 April 1997 


1007 


Keck I 


3910 - 8880 


10 


16 April 1997 


1012 


Lick 3 m 


3160 - 9870 


6 


11 May 1997 


1037 


Lick 3 m 


3120 - 9880 


6-11 


13 May 1999 


1769 


VLT 


3635 - 11680 


13 


12 June 1999 


1799 


VLT 


3908 - 9187 


5 



a Based on an assumed explosion date of 4 July 1994. 
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Table 2. 


Fluxes of narrow lines 






Line 


Rest 


Observed 




Flux a 








wavelength 


wavelength 22 June 6 Feb 


11 April 


12 June 


Notes 




(A) 


(A) 


1996 1997 


1997 


1999 




[re A1J 


7891.8 


7891.9 




1.2 


n i o 
U. i-Z 




7 




7610.3 




0.81 


U. 1 ( 




T A TTTl 

[Ar 111J 


7135.8 


7135.7 




0.26 


0.12 




TT „ T \ 

lie l yn) 


7065.3 


7065.5 






n ii 




fa. t n 
[Ar VJ 


7005.7 


7005.0 




0.21 




broad? 


J~1CK 


6562.8 


6563 8 




£ 2.1 


n vfi 
U. (o 




[le X] 


6374.5 


6375.5 


1.8 


1.3 






[U 1J 


6363.8 


6363 




0.47 






[U 1] 


6300.3 


6301.0 


3.0 


3.0 


n oo 




[re V11J 


6087.0 


6086 8 


3.7 


2.8 


n fin 
U.oU 




I'M TTl 

[N 11J 


5754.6 


5755.5 




2.2 




broad r 


[re VllJ 


5720.7 


5721.9 


3.0 


1.6 


n t 
0.7 




[re V1J 


5677.0 


5677.5 




0.50 






|"CV, WTTl i [TrU t 7Tl 

[be VllJ + [re V1J 


5276.4+5277.8 


5278.1 


6.5 


1.9 




bl. 


l"CV. WTl 

[re V1J 


5176.0 


5176.2 




3.3 


n ii 
U.ll 


Kl ,tt TTV, TT T7V-, "\7"TT 

bl. w. rC 11, rC Vll 


[he V11J 


5158.9 


5158.7 






n on 
U.2U 


Kl rr i „ tt tt„ t/t 
bl. w. he 11, he VI 


l"CV. WTl 
[re V1J 


5145.8 


5148.2 






n i o 

u. iy 


Kl ,tt TTV, TT TT^ "\7"TT 

bl. w. re ll, re Vll 


\r\ TTTl 

[U ill] 


5006.8 


5006.7 


44. 


26. 


o.O 




ro TTTl 
[U ill] 


4959.0 


4958.6 




8.5 


n so 
U.DZ 




T4 /3 


4861.3 


4861.6 




^3.6 


n ^1 




[IN e 1 V J 


4725.7 


4725.4 


1.8 


0.88 






[Ne IV] 


4714.2 - 4715.7 


4714.8 


0.43? 


0.41 






He II 


4686 


4686.0 


1.5 


0.97 


0.11 




[O III] 


4363.2 


4364.0 


6.5 8.3 


6.4 


0.53 




H 7 


4340.5 


4341.0 






0.27 




[Fe IV] 


4198.2 


4198.1 




0.07 






? 




4187.7 




0.05 






[Fe V] 


4180.6 


4179.8 




2.1 




bl. w. He I 4169 


[Fe V] 


4071.2 


4070.9 


3.0 


1.3 


0.11 




[Fe III]? 


4046.2 


4045.0 




0.04 




ident.? 


He I 


4009.3 


4009.0 




0.10 






[Ne III] 


3967.5 


3964.8 


3.4 






bl. w. He I 3964.1 


[Ne III] 


3868.8 


3868.1 


20. 9.5 


6.8 b 






[Fe V] + [Fe VII] 


3757.6+3758.9 


3757.7 


6.3 3.4 


9.0 b 




bl. 


[Fe VII] 


3586.3 


3585.3 


2.1 








[Ne V] 


3425.9 


3419.0 


15. 








[Ne III] 


3342.4 


3343.6 


3.8 









a Obscrvcd fluxes in units of 10 16 erg s 'cm 2 . 
b Actually measured in 11 May 1997 spectrum. 



-47- 



Table 3. Fluxes of intermediate lines 



Line 




Rest 
wavelength 

(A) 




Observed 
wavelength 

(A) 


22 June 
1996 


6 Feb 
1997 


Flux a 
11 April 
1997 


12 June 
1999 


Notes 


O I 




7771.9-7775.3 




7773.5 


1.4 




2.0 


1.1 




[O II] 




7319.9, 7330.2 




7326.5 


5.3 




4.6 


5.4 


bl. w. He I, Ca II 


[OI] 




6300.3, 6363.8 




6335.2 


4.5 




3.2 


4.4 




[O III] 




4959.0, 5006.8 




5006.7 


54. 




39. 


9.5 




Mg I] 




4571.1 




4574.1 


1.8 


1.3 


1.8 


1.6 


(+4583) 


[O III] 




4363.2 




4380.7 




5.3 


4.7 


1.2 


bl. w. H/3 


Mgl 




2852.1 




2862.3 




2.5 






bl. w. Mg II & Fc II 


Si III], C III] 


1882.7, 


1892.0, 1906.7, 


1909.6 


1899.3 




21. 






bl. 


N III] 




1746.8 - 1753.4 




1752.9 




2.1 








O III] 




1660.8, 1666.1 




1658.8 




9.9 








C IV 




1548.9, 1550.8 




1548.0 




7.6 








N IV] 




1483.3-1487.9 




1492.5 




2.2 






flux uncertain 


Si IV, O IV] 


1393.8, 


1402.8, 1397.2 - 


1407.4 


1406.3 




6.7 






bl. 


Si II, O I 


1304.4, 


1309.3, 1302.2 - 


1304.9 


1310.0 




2.1 






ident. uncertain 



a Observed fluxes are in units of 10 15 erg s 'cm 2 . 
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Table 4. Fluxes of broad lines 



Line 


Rest 


Observed 






Flux a 








wavelength 


wavelength 


22 June 


6 Feb 


11 April 


12 June 


Notes 




(A) 


(A) 


1996 


1997 


1997 


1999 b 





He I 


10830 


10827.6 








6.0 




Pe, He I,[S III]? 


9545.7, 9529.3, 9532.0 


9544.0 


3.0 






0.70 




PC 


9229 














Fe II 


9196.9-9218.2 










0.18 


bl. 


Fc II+Hc I 


9175.9, 9078.1, 9174.5 


9177.5 


7.0 






1.0 


bl. 


Fe II 


9122.9, 9132.4 


9133.8 


2.0 






0.39 


bl. 


Fc II 


9070.5, 9077.4 


9071 


0.54 






0.07 


bl. 


Pr? 


9014 










0.07 


bl. 


Ca II 


8662.1 


8668.2 






0.27 




bl. H I, Ca II 


H I 


8598.3 


8601.7 






0.07 




bl. 


H I, Ca II 


8545.4+8542.1 


8545.9 


0.20 




0.22 




bl. 


Ca II 


8498.0 


8496.9 






0.80 




bl. 


O I + Fe II 


8446.5 + 8451.0 


8448.1 


8.4 




3.5 


0.88 


bl. 


He I 


8285.4 


8287.3 


0.16 




0.12 


0.05 




Fc II 


8228.9 


8227.7 


0.72 




0.53 


0.11 




He I 


8155.8 


8156.3 






0.02 


0.02 




He I 


7971.6 


7970.7 


0.08 




0.22 


0.05 




? 




7892.0 






0.14 




narrow? 


? 




7610.4 






0.07 




narrow? 


He I 


7281.0 


7277.7 






0.46 




bl. w. [O II] 


Fc II 


7155.1, 7172.0 


7156.7 






0.21 




mult. 14 


He I 


7065.3 


7063.5 


2.0 




1.2 


0.40 




He I 


6678.2 


6675.6 


0.60 




0.33 






Ha 


6562.8 


6562.4 


330. 




160. 


30. 




[O I] 


6300.3, 6363.8 


6335.2 


4.5 




3.2 


4.4 




Fe II 


6148, 6149 


6146.3 






0.19 




mult. 74 


He I 


5875.7 


5876.4 


3.8 




2.5 


0.88 




[N II] 


5754.6 


5753.5 


0.39 




0.27 


0.11 




Fe II 


5534.9 


5533.2 






0.13 




mult. 55, bl. w. [O I] 5577 


Fc II 


5414.1, 5425.3 


5421.1 






0.18 




mult. 48+49 


Fc II 


5362.9 


5367.1 






0.22 




mult. 48 


Fc II 


5316.2 


5317.9 






0.54 




mult. 49 


Fc II 


5265, 5276.0 


5273.9 






0.10 




mult. 48+49 


Fc II 


5169.0 


5169.5 


2.6 




1.0 




mult. 42, bl. w. [Fc VI], [Fe VII] 


Fc II 


4629.3 


4627.6 






0.25 




mult. 37 


H/3 


4861.3 


4861.6 


8.8 




4.3 


1.2 




He II 


4685.7 


4686.0 


0.77 


0.24 


0.26 






He I 


4471.5 


4470.7 


1.1 


0.63 


0.68 


0.26 




[o in] 


4363.2 


4380.7 




5.4 


4.7 


1.2 


bl. w. H/3 


H7 


4340.5 


4341.0 




4.6 


3.4 


1.1 


bl. w. [O III] 


Fc II 


4233.2 


4238.9 






0.48 




mult. 27 


Fc II 


4178.9 


4180.3 






0.22 




mult. 28, bl. w. [Fc V] 4180 


He I 


4120.8 


4122.4 






? 






H<5 


4101.7 


4102.3 


1.7 




1.6 


0.78 




[Ne III], He 


3967.5, 3970.0 


3967.8 


2.2 




0.93 


0.30 


blend w. Fe V, He 


HC 


3889.0 














He I 


3888.6 




1.5 
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Table 4 — Continued 



Line 


Rest 


Observed 






Flux a 








wavelength 


wavelength 


22 June 


6 Feb 


11 April 


12 June 


Notes 




(A) 


(A) 


1996 


1997 


1997 


1999 


b 




3835.4 




1.2 


1.1 








[Nc V] 


3426 


3447 


1.8 










? 


? 

• i 


3075.4 




0.58 






uncert. id. 


? 


? 

■ > 


3001.1 




0.35 






uncert. id. 


O III 


? 

■ > 


2983.5 




0.36 






uncert. id. 


Fe II, [Mg V] 


2926.6 , 2928.3 


2925.8 




0.23 






mult. 60, ident. uncert. 


Mg II 


2795.5, 2802.7 


2793.3 




130. 








He I + Al II] 


2663.4+2669.2 


2666.3 




0.71 








Fc II 


2614.9-2631.0 


2617.7 




3.9 






mult. 1+171 


Fc II 


2585.9 -2593.7 


2592.0 




0.85 






mult. 1+64 


Fe II 


2562.5, 2563.5 


2561.3 




0.35 






mult. 64 


Fc II 


2506.4-2508.3 


2505.6 




2.0 








Lya 


1215 


1210.9 




300. 









a Observed fluxes are in units of 10 15 erg s 'cm 2 . 

b Includes some measurements from the 13 May 1999 VLT spectrum, especially in the ncar-IR. 



